function psi=computePsi(identityC, powersA)
psi = identityC * powersA;
end